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Abstract 



We consider the non-hermitian 2D Dirac Hamiltonian with (A): real ran- 
dom mass, imaginary scalar potential and imaginary gauge field potentials, 
and (B): arbitrary complex random potentials of all three kinds. In both 
cases this Hamiltonian gives rise to a derealization transition at zero energy 
with particle-hole symmetry in every realization of disorder. Case (A) is in 
addition time-reversal invariant, and can also be interpreted as the random- 
field XY Statistical Mechanics model in two dimensions. The supersymmet- 
ric approach to disorder averaging results in current-current perturbations 
of gl(N\N) super-current algebras. Special properties of the gl(N\N) al- 
gebra allow the exact computation of the /3eta-functions, and of the cor- 
relation functions of all currents. One of them is the Edwards- Anderson 
order parameter. The theory is 'nearly conformab and possesses a scale- 
invariant subsector which is not a current algebra. For N = 1, in ad- 
dition, we obtain an exact solution of all correlation functions. We also 
study the derealization transition of case (B) , with broken time reversal 
symmetry, in the Gade-Wegner (Random-Flux) universality class, using a 
GL(N\N; C)/U(N\N) sigma model, as well as its PSL(N\N) variant, and a 
corresponding generalized random XY model. For N = 1 the sigma model is 
shown to be identical to the current-current perturbation. For the dereal- 
ization transitions (case (A) and (B)) a density of states, diverging at zero 
energy, is found. 
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I. INTRODUCTION 



Quantum field theory in the presence of quenched disorder is an increasingly important 
subject in condensed matter physics. Unfortunately, even in two dimensions, fully solvable 
but non-trivial models are rather scarce. In this paper we consider, amongst others, what 
is perhaps the simplest model of disordered fermions that can be exactly solved but also 
exhibits some non-trivial behavior. We work in two spatial dimensions throughout. The 
model can be defined by the euclidean action^ 

St = J A-)^ L + ^ R {d z + A z )^ R + m{x) fy L +m\x) ^ R ) (1.1) 

where d 2 x = dxdy, z = (x + iy)/2, z = (x — iy)/2, and A z = A x — iA y , A z = A x + %A y . The 
subscripts L, R on the fermion fields denote left versus right movers in the conformal theory 
with m = m* = 0. Letting m = V + iM, m* = V — iM, we take V and M to have gaussian 
probability distributions of equalQ strength g and zero mean, and similarly for A x , A y with 
strength gx- 

P[m,m*] = exp f — f m*(x)m(x)\ (1-2) 

V 9 J 2ix I 



P[A Z ,A-] = exp (-— [ ^ A z (x)A z (x 



(The couplings g,gA represent twice the variance.) Disorder averaging corresponds to a 
gaussian functional integral: 

JO) = J Dm*DmDA z DA- P[m, m*}P[A z , A-} (O) (1.3) 

In the context of two-dimensional statistical mechanics, the above model may be thought 
of as a random field XY model or equivalently the random field sine-Gordon model M, at 
the free-fermion point. The usual sine-Gordon model has m constant, which now acquires 
a random phase, preventing the theory from becoming massive. Random field XY models 
were first studied by Cardy and Ostlund [0] using replicas ||. The action in Eq.( |l.l| ) can 
be related (see below) to a 2D Dirac hamiltonian H 2 subject to a real random mass M, 
a random imaginary potential — W with strength gu = 9v = 9, and also in a random 
imaginary gauge-field 

H 2 = {-id x + %A x )a x + (-id y + iA y )a y + M(x, y)a 3 - iV(x, y) (1.4) 

where Oi are the standard Pauli matrices. This hamiltonian is non-hermitian. This theory 
should be compared with the models described in Ref.'s 0], || which were argued to have a 



1 partition function Z = J T>[ifj^,ip] exp(— S) 

2 The difference of the variances of Gaussian random variables V and M renormalizes to zero. 
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transition in the universality class of the quantum Hall plateau transition; the latter models 
have a real random mass and a real scalar potential with strength gu ^ gy, and a real 
random gauge field, the corresponding 2D Dirac hamiltonian being hermitianP]. The same 
model as in Eq. ( |1.4|) , also describes^ the continuum limit of electrons hopping on a square 
lattice with tc flux ||, exhibiting a derealization transition at zero energy, first discussed 
by Hatsugai, Wen and KohmotofHWK] in Ref. [§]. This involves a doubling of the number 
of degrees of freedom( see section II and Appendix B), leading to a random hermitian 4- 
component hamiltonian. Applications of our results to this derealization transition will be 
described in section V . 

We perform the disorder averaging using the supersymmetric method. This leads to a 
left-right current-current perturbation of the gl(N\N) conformal supercurrent algebra for the 
iV-species version, used to compute Nth moment disorder averages. Conformal field theory 
methods were previously used by Bernard for the case where guipi gv) — [||, which 
leads to the osp(2N\2N) algebra. Also, the case of gy = gu = 0, but gA ^ 0, where the 
conformal symmetry is unbroken and the model is equivalent to free bosons, was studied in 
0,0. Though our model has a non-zero renormalization group /3eta-function, some special 
features of the gl(N\N) algebra, in particular the nilpotency of the fermionic generators of 
gl(N\N) and the existence of two quadratic casimirs, lead to a complete solution for all 
correlation functions in the case N = 1, and to the solution of the correlation functions of 
all Noether currents in the case of general N. The only non-vanishing /3etafunction is that 
for gA which we compute exactly. It only depends on g. Integration of the renormalization 
group equations allows us to compute for instance the density of states of the derealization 
problem of ||. 

Our results reveal some interesting features of the current-current perturbations of Lie 
superalgebras in comparison with the analogous perturbation of the ordinary bosonic Lie 
algebras[|. Our models are integrable quantum field theories for the usual reasons (Lax pair, 
etc.), and are not conformal. However the structure of the theory allows a solution which 
doesn't rely on exact S-matrices, form-factors, etc. Regarding the difficulty of solution they 
lie somewhere between the conformal field theories and the massive integrable theories. 

Derealization transitions of 'sublattice' (or 'random flux') models, studied first by 



Hikami et al. [11], and by Gade and Wegner [E| perturbatively using replica sigma models, 



have been known for some time to exhibit R.G. beta functions of the kind we obtained for 
the current-current perturbations of gl(N\N) supercurrent algebras. These derealization 



3 We point out that for this Quantum Hall setting, based on the /3eta functions computed in ||, 
the line gM + gv = which is the subject of this paper, is a line of fixed points in the ultraviolet. 
One can check numerically that starting from gM, gv, 9 A all positive, one can actually flow to the 
line gM + <?y = in the UV. 

4 We were first made aware of the relationship between the model of Hatsugai et al. and the 
random XY model by M.P.A. Fisher in the context of the replica field theory. 

5 The current-current perturbation of su{2) at level 1 gives the Gross-Neveu model, which is a 



massive theory with factorizable S-matrix. See for instance [ 10 1 . 
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transitions possess a special particle-hole symmetry, known to prevent localization due to 
disorder even in one dimension This is not an accident: In the last section of this 

paper we exhibit a random Dirac Hamiltonian giving rise to these theories (in 2D). From 
this Dirac Hamiltonian we derive a SUSY sigma model, which we find to live on a target 
manifold GL(N\N; C)/U(N\N). We solve this sigma model exactly for the case N — 1, in 
the presence of a topological Wess-Zumino-Witten(WZW) term with arbitrary real coupling 
k. We find that the GL(1|1; C)/U(l\l) sigma model with WZW term is the same theory as 
the current- current perturbation of gZ(l|l) supercurrent algebra. One expects that such an 
equivalence should extent to the case of general N. 

An interesting and special feature of our current-current perturbations is that (for any 
value of N) the non-vanishing /5etafunction is associated only with two of the generators of 
the gl(N\N) algebra. These can be separated out and a scale invariant fixed point theory 
results for any value of the disorder strength g, i.e. a line of fixed points (for any value of N). 
After we obtained the current algebra results in section IV, the papers flM |T3 appeared 



on PSL(N\N) sigma models, exhibiting a line of fixed points as the sigma model coupling 
constant is varied. We have added a paragraph to section IV C explaining the connection 
of our results with these papers. Our bosonized formulation may provide a useful way of 
studying the PSL(N\N) sigma models. Also, a proposal for a relationship between the PSL 



sigma models and the integer quantum Hall plateau transition has appeared recently [16 



The outline of the paper is as follows. In section II we give a general symmetry-based 
overview of the models we discuss in this paper. Sections III and IV address the random 
field XY model. In section III we formulate the supersymmetric method and extend it to the 
iV-species version. In section IV A we compute the /5eta-function for the random field XY 
model to all orders, and the exact correlation functions using current algebra techniques. In 
section IV B some of these results are extended to the A-species case, possessing gl(N\N) 
symmetry. The scale invariant subsector of this theory is discussed in more detail in Section 
IV C. A generalized random field XY model, whose associated Dirac hamiltonian lacks 
time-reversal symmetry is discussed in Section IV D. (This Dirac Hamiltonian underlies 
the Gade-Wegner 'sublattice' models discussed later in section VI.) In section V we turn to 
derealization transition of Hatsugai et al. || and compute the density of states as a function 
of energy. In section VI A we discuss the corresponding derealization transition with broken 
time reversal symmetry, which is in the Gade-Wegner universality class |12|] and find that it 
is described by a GL(1\1; C)/U(l\l) sigma model, which we solve exactly. In section VI C 
we relate the sigma model to the current-current perturbation discussed earlier (in section 
IV A). In section VI C we conclude with a summary of the A r -species generalization, the 
GL(N\N;C)/U(N\N) sigma model, and its scale invariant PSL(N\N) variant. A number 
of technical details are delegated to three Appendices. In Appendix A we relate the N — 1 
species theory discussed in section IV A, to the random field XY model. Appendix B gives 
the steps leading from the general non-hermitian Dirac hamiltonian and its hermitianization 
to a corresponding conformal field theory action. Finally, Appendix C contains details of 
the 1-loop RG equations of Section IV D. 
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II. THE MODELS AND THEIR SYMMETRIES: AN OVERVIEW 



The quantum mechanics of non-interacting particles (electrons) subject to random po- 
tentials and their possible universality classes of localization/delocalization transitions are 
conventionally classified according to those symmetries which are present for a fixed realiza- 
tion of disorder. In this section we describe the class of models discussed in this paper, all 
of which exhibit particle-hole symmetry. 

A. Time Reversal Symmetric Models 

The fundamental model that we study in this paper can be formulated as the random 
field XY model |2|. This is a 2D statistical mechanics model. At the free fermion point of the 
XY model, however, the XY spin operator can be represented, using abelian bosonization, 
by a bilinear ipL^R °f Dirac fermions governed by the action of Eq. (|1.1| ). Alternatively, this 
action can also be related to the 2-component random Dirac Hamiltonian H 2 of Eq.( |1.4j ) in 
two spatial dimensions, containing a random real Dirac mass M, a random imaginary scalar 
potential V and purely imaginary gauge potential A term. (Details of this connection are 
lain out in Appendix A) . The latter being a non- interacting system, this allows us to use the 
supersymmetry (SUSY) method for disorder averaging. This is the method we use in this 
paper to study the random XY models. An apparent difficulty in applying this method is 
the non-hermiticity |HJ of this quantum mechanical hamiltonian if 2, which seems to prevent 



the functional integral over the bosonic fields, appearing in the SUSY approach to disorder, 
from converging. This problem is however easily remedied by considering an associated 
hermitian 4-component hamiltonian (Appendix A) 



\Hl 

As discussed by HWK || , the hamiltonian H± itself describes the quantum mechanics (in the 
continuum limit) of fermions hopping, with real hopping amplitudes, on a 2D square lattice 



with 7r-flux through each plaquette (see also [18]). This theory possesses a derealization 



transition at zero energy, the characteristic symmetries of which are particle-hole and time- 
reversal symmetries^. Hence there are two apparently unrelated problems linked to this 
4-component hamiltonian, the random field XY model, and the lattice fermion hopping 
model. In the latter model, a small imaginary part has to be added to define the quantum 
mechanical Green's functions. This is not necessary for the interpretation as an random XY 
statistical mechanics model. As discussed in more detail below (and in Appendix A), the 
theory describing averaged quantities of the random XY model has gl(l\l) global SUSY, and 
gl(N\N) SUSY when higher moment averages (e.g. of 'Edwards- Anderson order parameter' 
type) are considered. The corresponding derealization transition ||, on the other hand, 
can be shown to be a GL(2\2; R)/Osp(2\2) sigma model. 



6 For the action of these symmetries in the continuum theory with hamiltonian H4, see Eq. 

dHJH). 
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B. Models with Broken Time Reversal Symmetry 



In the random field XY model discussed above, the complex XY spin, represented by 
iI>li/jr is coupled to a random ordering field, which is a random phase m(x) = me"' with 
•§ real. The corresponding random 2D Dirac Hamiltonian H 4 ( of Eq. (|1.4j , |2.1|) ) is time 
reversal invariant. Time reversal symmetry breaking may for example be introduced by 
allowing, in the language of the Dirac Hamiltonian H 2 , in addition to the purely imaginary 
gauge potential discussed above, also a real gauge potential^ The addition of a real random 
gauge potential leads to the generation of both, real and imaginary parts of all three types 
of potentials, and the theory flows off to strong coupling^]. It is then more useful to employ a 
description in terms of a non-linear sigma model (section VI). The corresponding quantum 
mechanical hamiltonian if 4, possessing particle- hole symmetry but not time- reversal sym- 
metry, is now, by universality, in the class of the Gade-Wegner sublattice models [0, |19 



describing fermions hopping on the square lattice with complex hopping amplitudes. The 
corresponding random XY model has GL(1|1;C) global SUSY. In section VI we describe 
the corresponding Gade-Wegner type derealization transition in terms of a non-linear sigma 
model on a manifold which we denote by GL(1|1; C)/U(l\l) (allowing also for a topological 
Wess-Zumino term.) 

It is noteworthy to re-emphasize that both models, with or without time reversal sym- 
metry, possess particle-hole symmetry at zero energy, a feature which can circumvent local- 



ization of all states due to disorder, at that energy, even in one dimension W% 



III. SUPERSYMMETRIC DISORDER AVERAGING 

In this section and the following sections IV A, B, C, we discuss the random field XY 
model, described by the hamiltonian of Eq.( |1.4| ), which we can solve exactly. 

To implement the supersymmetric method for disorder averaging ||20|| , we augment the 
theory with bosonic ghosts /3lr,/3l,r, coupled to the disorder in the same way as the 
fermions, with action S$ = Sq(i() — > 0). The partition function of the fermions plus ghosts 
is now independent of the disordered potentials^ so the order of functional integration over 
the 'matter-ghosts' and disorder can be interchanged. Integrating over disorder, this yields 
an effective action for the fermions and ghosts: 

ScS = I 27 (W^ £ + ^ 9 ^ R + + PR 9 *?*) + 9°m + 9aO a ) (3.1) 



7 Properties of the most general Dirac hamiltonian where random Dirac mass, scalar and gauge 
potentials have both real and imaginary parts, as well as the relationship with a corresponding 
SUSY field theory, are summarized in Appendix B. 



See Eq.( 4.84| ) of section IV D, with the notations of Appendix B 



9 There is a subtlety arising from the convergence of the bosonic functional integral. A careful 
discussion of this issue, done in Appendix A, results in the effective action described below. 
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where 

o m = - (i>y L + &pl) tylfa + pIpr) (3.2) 

o A = - (^l + PlP L ) (i>y R + P R p R ) 

In the sequel we will need to introduce a multi-species version of the above model. 
Introducing a species index ipLa, a = 1, •-, N, we define the iV-species model as 

f d 2 x N 

Si = / -FT E VPlJPs + A ~)^La + Wlaidz + A)^Ra + m(x)^L^La + m*(x)^L^i?a (3.3) 
J 271 a=l 

Introducing bosonic ghosts and integrating over disorder leads to the effective action: 

Sis } = J ^ (E (^L^La + V>L^ila + /3L^ia + /^Aftta) + gO m + 9a O^J (3.4) 

where now 

O m = -J2 {lhl>La + PkPLa) (^Rb + P{ b P Rb ) (3.5) 
a, ft 

Oa = -Y1 {4J>l* + PLPlo) {il>y>m + P\u>Pm) (3-6) 

a, ft 



IV. RANDOM XY MODELS VIA GL(N\N) CURRENT ALGEBRAS 

A. Random Field XY Model: N = 1 Species 

In this section we consider the N = 1 model in detail. The bosonic ghosts /3l (Pr) have 
Lorentz spin 1/2 (—1/2). The conformal field theory for such a bosonic first order action 
was treated in ref. |21|. The Virasoro central charge for one copy of the ft', ft system is 
c = —1, so that the total central charge of the matter plus ghosts is zero. This had to be 
the case since the partition function is one. 



In the effective action eq. ( 371]) for N — 1 species the operators O rn and Oa are bilinears 
in the dimension 1 currents: 

_G+ = ftkl> L , _G-=Pl^ l , _J=:Ml--, J=--PlP{: (4.1) 
G + = P R ij R , GL = /3 fl Vt, J=:Vfl4:, J'=-PrPr- (4-2) 

(We denote all left-right currents as (Jl, Jr) = (J, J).) Namely, 

O m = JJ- J'7' - GJG + + (4.3) 
O a = -(7 -J)(f -J) 



In the conformal limit g — g A — 0, the currents generate a g/(l|l) current algebra. The 
currents J, J' are bosonic Z7" (1) currents for the fermions and ghosts respectively, and G± 
are their fermionic superpartners. From the operator product expansions (OPE) 
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M*WlW = ^l(z)M^) ~ j^, PIWPlW = -p L (z)Pl(w) ~ - j— (4.4) 
one obtains the OPE's: 

J(z)J(w) ~ 7—^5, J'(z) J'H ~ -7-^5 (4-5) 

J(z)G±H ~ ± G±(w), J'(z)G±(w) ~ ±^^G±(w) 

2 — U7 2 — W 

g_(*)g + m ~ 7 -^— - + — (jh - J'M) 

(z — u>) 2 z — w 

Let j,j',g± denote the zero modes of the currents, defined as j = § j^J(z), etc. These 
zero modes are generators for the global g>i(l|l)algebraR: 

V, g±\ = b", a±\ = ±g±, {g+, g-}= j- f, g± = o (4.6) 

The algebra g/(l|l)has two quadratic Casimirs: 

C = f - f 2 - g-g + + g + g- } C = (j - j'f (4.7) 

The operators O m and O a are left-right current- current perturbations with precisely the 
structure of the Casimirs C,C respectively. This implies that the gl(l\l)i, ® gl (1\1)r sym- 
metry of the conformal field theory is broken to a diagonal g/(l|l)symmetry in S^g. 

The model defined by S e s is a perturbation of a c = conformal field theory, and can 
be studied in the framework of Zamolodchikov [22]. The g/(l|l) conformal field theory was 



studied in [J23|. Current-current perturbations are generically integrable in the case of the 



bosonic Lie algebras, which means there are an infinite number of conserved currents. It is 
likely that this feature also holds for the current-current perturbations of the superalgebras. 
However we do not pursue this here since, as we now show, the model can be solved by 
straightforward manipulations of the functional integral. 

The manner in which conformal symmetry is broken in S e g is contained in the ^eta- 
functions. The OPE's of the perturbing operators are the following: 

O m {z,z)O m {0) ~ ^O A (0), O A {z,z)O m {0) ~0, O A (z,z)O A (0) ~0 (4.8) 
Note that because of the supersymmetry, the l/(zz) 2 terms vanish. One obtains to lowest 



order [g4 



Pa(h) = ^ = 29 2 + .~ (4-9) 

(3 g = 



10 For the complete set of modes the current algebra (4.6) defines the affine Lie superalgebra 
gl{l\ 1). 
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where I is the log of the rescaling factor. The ultraviolet limit corresponds to I — > — oo 
where — > — oo, whereas the infrared limit is I — > oo with g^ — * +00. This shows that if 
qa disorder were not included from the start it would be generated under renormalization. 
Since the operator O m is never generated in the OPE at higher orders, it follows that j3 g = 
to all orders. 

It will be useful to bosonize the theory. The U(l) fermion current is bosonized with a 
scalar field <fi: 

J = id z <p, J = -idz4> (4-10) 

In the conformal limit, (f)(z,z) = (pi(z) + <Pr{z) an d 

"0L = exp(zy9 L ), ip{ = exp(-i(p L ), ip R = exp(-i(p R ), ip R = exp(i(p R ) (4.11) 

Bosonization of the ft — (3 system was described in The J' currents are expressed 
in terms of another scalar field (f)'\ 

J' = id K <j>', 1' = -idsfi (4.12) 

The scalar field <$' has the opposite sign kinetic term compared to <ft, which implies: 

<p(z, z)<f>(0) ~ - log zz, <p'{z, z)<p'(0) ~ + log zz (4. 13) 

Since the 0' field still contributes c = 1, in order to obtain the c = — 1 of the ft — (3 system 
one needs an additional fermionic 77 — £ system with c = —2. The result is 

ft L = e -^d z £ L , p L = e l ^ VL , ft R = e^d^ R , (3 R = e" 4 ^ (4.14) 

where as before 0' = ip' L + ip' R . 

The r) — £ system has the first order action 

Sn-i = / Vl^ l + rj R d z £ R (4.15) 
J Ztt 

In order to have a more symmetric effective action it is convenient to formulate the 77 — £ 
system in terms of a complex dimension zero fermion X i with the action [25] 

S x = J^ (d zX %X + ^X ] d zX ) (4.16) 

The equation of motion d z &zX = implies that X separates into left and right moving parts: 
X — Xl{ z ) + Xr(z)- The 77 — £ fields are related to X as follows: 

d z £ L = id zX , dz£n = idzX (4.17) 
Vl = id z x\ Vr = idzX ] 

In order to absorb the — J J + J'j' terms in O m into the kinetic energy, we define 

*± = \{*±f), (418) 
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Putting all of this together one obtains the effective action 
d 2 x r 



ScS = I 47 l 9 '^ 9 ^- ~ 2 C 2 9Ad z (t>-^~ + (l + 2ge^-) d zX %X (4-19) 
+ (l + 2ge«*-)d sX %X 



In order to compute correlation functions, we treat the effective theory in a number of 
steps. First note that 



S, 



efi 



, 9 a] = S, 



cff 



2CW- 



9a = 



(4.20) 



Let $(0+, 0_, x) De a composite field or product of such fields at different space-time points. 
Making a change of variables + — > (f> + + 2( 2 gA<fi- in the functional integral, one relates 
correlation functions with g& ^ to those in the qa = theory: 



(4.21) 



Henceforth S e s denotes the effective action (|4.20|) with qa = 0. 

The next step recognizes that <p + acts as a Lagrange multiplier. Now let $(0-, x) denote 
a field or product of fields that is independent of <p + . Introducing a source p for + then 
functionally integrating over + gives a functional 5 function: one has 

e -/€r^+$(0_, x )\ = / J D[0_, X ]5(^0_-p)e-^^=°]$(0_,x) (4.22) 

/ g A =0 J 

= J D[x\x}e- S ^ +=0 '^ = M<I>- = M 



where p is the solution to 



dz&zp = P 



(4.23) 



(We have set the partition function to 1.) 

Insertions of the operator d z dz<f) + in a correlation function can be obtained from func- 
tional derivatives with respect to p. For instance, 



One needs 



-47T 



5p(z, z] 



J D X e-^^=°^-=^$(0_=p,x) 



Jp=0 



4vr— S'efff 
5p 



0. 



P] 



p=0 



X 



and 



S eff [0 + = 0_=O] = (l + 2<?)^ 



We define a rescaled x field, 



(4.24) 



(4.25) 



(4.26) 



(4.27) 
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such that 



S~ = (1 + 2g)S x = J ^ (Ux^X + <tfd. 



(4.28) 



The above formulas imply that all correlation functions are reduced to a free-field corre- 
lators in the rescaled x system defined by S~. For instance: 



9A=0 



- Att 



s*(<f>-,x) 



(4.29) 
(4.30) 



<t>-=oi 



All correlation functions ( )j on the RHS of the above equation are computed with the free 
action \ (|4.16| ). In particular 



X j (z,z)x(0)) = ~ (x(z,z)xHty) = -log(^) 



(4.31) 



The second term in eq. (|4.29|) generally involves 5 functions. More generally the insertion 
of d z dj(p + at different space-time points is equivalent to the insertion of the operator on the 
RHS of eq. ( |4.25| ) at the same points in the free theory x U P to 5-functions: 

((d z «3,0 + r$(0_ x )V = (2igC) n ((d z x%X-^x i d z xT^ = 0, x ))~ + 5 - functions 



(4.32) 

The above formulas are sufficient to compute all correlation functions. First, for corre- 
lations not involving + , the field 0_ may be set to zero: 



(*(0-,x)> 



9A=0 



Thus, 



Equation ( 4.29|) implies: 



Thus to all orders in g: 



(z,*)0_(O)) = O 



d z d^((f) + (z,z)(f>^) gA=0 = -Atx5 2 (x) 



9A=0 



— 2 log(,z,2:) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



Correlation functions involving multiple 0+ insertions are more interesting. From (|4.32|) 
one finds 



(d z dz<f) + (z, z)d w &uj(j)+(w,w)} = -8£f 2 C 6 - 
Integrating this, one finds a log 2 contribution: 



1 



z — w) 2 (z — w) 2 



+ 5 — functions 



(4.37) 
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(<M^)<MO)> 9A=0 = -VC 6 log 2 (^/a 2 ) (4.38) 

where a is an ultraviolet cutoff scale. 

Finally let us compute correlation functions of the ' matter' U(l) field <fi = C(0+ + 0-)/2 
with g A 7^ 0. Using eq. (|4.21| ) one obtains 

(<f>(z,z)4>(0)) = "C 2 (l + 2( 2 g A )\og(zz/a 2 ) - g 2 ( 8 log 2 (zz/ a 2 ) (4.39) 

Similarly, 

(<P'(z, z)<f>'(0)) = C 2 (l - 2C V) log(^/a 2 ) - /C 8 log 2 (^/a 2 ) (4.40) 

The /3eta- function /3 A (g) to all orders in g can be computed from (|4.39|) along with the 
renormalization group equation 

/ ° +PA(g)^-)(<P(r)m) = (4-41) 



dhga dg A 



One finds 



(3a(9) = 2g 2 ( 4 = {1 2 + g 2gr h = (4-42) 

Note that this agrees with the lowest order result ( |4.10| ). 

We consider now correlation functions of exponentials of cf) and <$' . The 'free-field prop- 
erty' eq. ( |4.32| ) implies that correlation functions of exponentials can be computed using 
Wick's theorem in the free theory S x . One has 

/jje^^A = Hexp (-a iaj (0(^)0(^-)>) (4.43) 

\ i I i<j 



n \zi - Z3 \ 2 «^ + ^ gA) exp ^ a . a . g 2 C s log2 {z , _ (444) 



i<j 

From this one deduces how the disorder affects the dimension of the exponential fields: 

dim ( e **) = a 2 C 2 (l + 2( 2 g A ), dim (W) = -a 2 C 2 (l - 2( 2 g A ) (4.45) 

B. Random Field XY Model: iV-Species 

We now extend some of the results of the last section to iV-species. The effective action 
eq. ( |3.4|) has a gl(N\N) Lie superalgebra symmetry. To see this, define the left-moving 
currents 

J ab = ML J U = Pirfh (4-46) 
G ab = pU La , Gl = Ml 

and similarly for the right-moving currents J a b = ipmWRai e ^ c - 111 the conformal field theory 
with g = g A = 0, one has a gl(N\N) current algebra with the operator product expansions 
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Jab{z)Jcd{w 

J'ab(z)J'cd( W 
Jab{z)G cd {w 
Jab{.z)G ] cd {w 

Gab{z)G\ d (w 



k 5 ad 5 bc 
(z — w) 2 
k 5„a5, 



z — w 



ad^bc 



(z — w) 2 z — w 



(5adJbc(w) - hcJad{w)) 

{8 ad J' hc {w) - 5 bc J' ad {w)) 



5 n 



z — w 

he 



G bd (w) 



J'ab(z)Gcd(w) 



5, 



Id 



Z — W 



■G ca (w) 



k 5 ac S bd 



5, 



ad 



1 



(z — w)' 



z — w 



z — w 

{8acJbd( W ) ~ SbdJcaiw)) 



GlW 



(4.47) 
(4.48) 
(4.49) 
(4.50) 
(4.51) 



where the level of the current algebra is k — 1. 

The perturbing fields are current-current perturbations which preserve the global diag- 
onal (left-right) gl(N\N) symmetry: 



a 



EG* 

a, b 



baJ ab 



r i' 

°ba° ab 



^ab^ab - ^ba U 



ba 



(4.52) 



O 



A 



E (J'ac 
a,b 



Jn 



J bb 



J 



bb 



The structure of the perturbative expansion in g, about the conformal gl(N\N) current 
algebra closely parallels the 1-species case. In particular, using the OPE's ( |4.51| ), one can 
verify that the OPE's ( |4.8|) still hold. In particular, the iV dependence cancels. Since the 
ultraviolet divergences arise from the singular terms in the OPE eq. ( [4.8[ ), this implies 
that the /3eta- functions are independent of N and are the same as in eq. ( |4.10|) . Thus the 
N-species theory has the same /?eta function as the 1-species theory. 

Since the perturbative expansion can be derived from ( |4.8| ), we expect that correlation 
functions in the gl(N\N) theory have the same structure as we found for gZ(l|l). However, 
since the functional integral manipulations we used to solve the bosonized gl(l\l) theory 
are not straightforwardly extended to the gl(N\N) case, here we obtain expressions for the 
correlation functions by using properties of the current algebra. This can be viewed as an 
alternative method to the functional integral methods described in section IV A, and in fact 
explain the results we obtained there in a more covariant fashion. 

Let J a denote an arbitrary left current, J a G {J a b, J'abiG a b,G ab }. The perturbing opera- 
tors O m , Oa are gl(N\N) invariant and can be expressed as 



C a0 J n J 



a<J 0i 



a 



-C af3 J n J 



(4.53) 



where C a/3 and C af3 define two independent Casimirs. In the abelian sub-vector space 
spanned by (J u , J 22 , Jnn, A 



., J' NN ), one has 



{C al3 } 



In 






-In 



{C af3 } 



{1} -{1} 
-{1} {1} 



(4.54) 



where Iat is the N x N identity matrix, and {1} is the N x N matrix with l's in every 
entry. Since C is invertible (whereas C is not), C defines a metric to be utilized to raise 
and lower indices. We will need C a p = C aa rCppiC a 13 , where C a p = C alS . In the subspace 
defined above, 
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{C a p} -- 

The OPE's ( 4.51| ) can be expressed as 

JJz)Jr(w) ~ 



{1} {1} 
{1} {1} 



(4.55) 



a 



a/3 



+ 



(z — w) 2 z — w 

where /Zg are structure constants. 

The U(l) currents can be bosonized as in ( |4.10| , |4.12 ): 



-fa/3 J i( w ) 



(4.56) 



Jn 



Jn 



-id^f) a , J' aa = id z <f)' a , J\ 



aa 



(4.57) 



By bosonizing the fields as in (4.11, 1.14| ), one can verify that the shift property ( [4.2(J ) 
continues to hold: 



S e S [4>a, 4>' a , 9a] = S cS (j) a - (V^, 4>' a ~ C^S, fiU = 



where 



a 

and ( is defined in ( |4.18 ). This implies 



9A=0 



(4.58) 



(4.59) 



(4.60) 



where O is any operator. 

The field £ has very simple correlation functions because it doesn't appear in any ex- 
ponential interactions. Collecting the kinetic terms for the scalar fields coming from the 
conformal field theory and from the terms g(J aa Jaa — J'aaJaa) 



1 r d?x a 

^kinetic = J — d z<K&z<K 



where 



(4.61) 



(4.62) 



Since £ doesn't couple to any of the remaining interaction terms, one has the exact 2-point 
function 



#(*,*)£(o) 



-2C 2 log zz 



(4.63) 



The above properties, together with current conservation provide strong constraints on 
the 2-point functions of the currents. Euclidean rotational invariance, and the fact that 
J a , J a are dimension 1 currents constrains the (J J) 2-point function to be of the form: 



J a (z, z)Jp(0)) = — (C af3 F(g, zz/2 ) + C af sF(g, zz/j,' 

' 9A U 2,2, ^ 



(4.64) 
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where we have used the unbroken gl(N\N) symmetry. F,F are scaling functions of g and 
the dimensionless combination z~zfi 2 , where \x is the energy scale that enters through the 
renormalization group. 

The field £ has the covariant description: 

id z Y, = & a Jp, -id^ = &J^ Vq (4.65) 

From ( |4.63|) we see that away from z~z = 0, one must have C^,(J a Jp) = 0. Since C^,C a p ^ 0, 
and C^,C a p = 0, one deduces that F — 0. 

The property (|4.63|) , together with (|4.60|) , leads to g& independence (up to ^-functions) 
of the correlation function ( |4.64j ). Since the dependence on the scale [i must be compatible 
with the /3eta-f unction for gA, the qa independence implies that F is only a function of the 
dimensionless coupling g. Thus 

(j a (z,z)Jf,{0)) gA=o = ^f al3 F{g) (4.66) 

Current conservation, d% J a + d z J a = 0, requires 

dz(J a (z,z)J p (p)) = -d z (j a (z,z)Jp(0)) (4.67) 

Integrating this, 

(J Q (z, z)M0)) gA=0 = ^ (C 2 C aP + C aP F{g) \ogzz) (4.68) 

The coefficient ( 2 of the C a p term was derived from ( f4.63|) . 

To fix the one unknown function F(g), we use current algebra equations of motion. The 
equations of motion to first order in g can be obtained in conformal perturbation theory 

r dw 

dzJ a (z,z) = g j —J a (w,z)O m (z,z) (4.69) 

Taking into account the £-rescaling — > <p/( that gives the kinetic term the standard 
normalization, and using the OPE's ( 4.56|) , one obtains 

dzJ* = 9( 2 ft~<Jf3J, = ~9ja (4.70) 
Using this in the correlation function (dzJ a d z J p) , along with ( |4.64j ), 

J^C^F = gVft'y a / ((J^)(z,z)(J a Jy)(0)) gA=Q (4.71) 

The non-zero contribution to F comes from 

((Jp,J,)(z,z)(Jj Y )(0)) ~^Cp, a C w (4.72) 

x ' 9 A— <J Z Z 

where we have used ( [4.681 ). Now, in terms of the structure constants, the OPE 
O m (z,z)O m (0) 20 A /zz implies 
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(4.73) 



The last three equations thus give 



2^8 



F(g) = 2gX 



(4.74) 



Let us summarize the above results by integrating the current correlations to find the 
2-point functions of the scalar fields: 



0+(z,Z)#(O)) ~ -4g 2 ( 8 log 2 zz, V a,b 

I 9.4=0 

<j>t(z,z)fc(0)) ~-2( 2 5 ab logzz 

I 9.4=0 



' 9A=0 



Finally, using the shift property ( |4.60|) , 

(0«Cz,s)&(O)) = -C 2 (<U + 2C 2 g A ) \ogzz-g 2 ( 8 \og 2 zz 



(4.75) 
(4.76) 
(4.77) 

(4.78) 



U*,*Wb(P)) = C 2 (<U ~ 2( 2 9A ) log zz -g 2 ( 8 log 2 zz 



(4.79) 



C. Nearly Conformal Structure 

We have seen in the previous section that the N-species theory (describing N copies of 
the random field XY model) is nearly conformal. This can be seen more clearly by sepa- 
rating out the non-conformal pieces (coupling holomorphic and anti-holomorphic coordinate 
dependences) by rewriting Eq.'s( fO)6| , |4.68| ) as follows. Instead of the basis of currents as in 
Eq.( [1.47|) , we extract U(l) x U(l) from the bosonic U(N) x U(N) subgroup. Explicitly, 
consider traceless generator matrices T A of su(N), (A = 1, N 2 — 1) 

J A = (TA)abJab, J' A = (T A ) ab J' ab 

and 

J = ^ t J aai J = ^ t J aa 

a a 

and similar for the fermionic currents (which we denote by J A , A = 2N 2 + 1, AN 2 ). 
Forming the combinations 

J_ = J - J', J+ = J + J', 

and noting that all non-conformal terms in the current-current correlators are multiplied by 
the invariant Cab whose only non- vanishing matrix elements are A, B = ± in this basis, we 
see that the two point correlation functions of all the remaining currents are 

< J A (z,z)J B (0) >= < J A {z,z)J±{0) >= (A,B^ ±) 
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Note that we see from Eg. ( |4.72| ) that these correlators are invariant under independent 
action of the symmetry group on the left- and right- chiral components of the currents. 

The ability to separate off the J± currents (and the corresponding fields, 0±) can be seen 
from the hamiltonian, expressed as 

H = H + H I 

where H Q is the non-interacting g/(iV|iV)-Sugawara Hamiltonian, and Hj is the current- 
current interaction Recalling that Cab = C AB there are no terms in the hamiltonian H 
which couple the generators J±, J± to the remaining ones, which allows to consistently set 
J± = J± = 0, yielding a scale invariant theory. 

Setting J± = J± = effectively sets g& = 0. Since the non-zero /3eta-function is P 9A (g), 
what is left is a conformal field theory since f3 g = 0. In terms of groups, setting J± = is 
equivalent to dividing by the £7(1) ® £7(1) subgroup: GL(N\N)/U(1) ® £7(1) = PSL(N\N). 
Since the two £7(1) bosons give c = 2, and the total central charge for GL(N\N) is zero, 
what is left is a c = —2 conformal field theory. The conformal PSL(N\N) sigma model was 
recently studied in |L4] Jl5| (see also section VI below). 



In the N = 1-species case the emergence of conformal symmetry is easily understood: 
Bosonizing the non-interacting theory, there are only two generators G + = idx, G_ = idx 
left once the currents J± and thus the fields <f>± have been eliminated. The PSL(1\1) sigma 
model is a free theory. 

One proceeds similarly for the N > 1 cases. The N = 2 species case for example describes 
the second moments of observables in the random XY model. An interesting observable is the 
'Edwards- Anderson order parameter', which is a bilinear of one of the 'conformal' currents: 

912 = (*pLl^Rl)(^L2^R 2 ) = Jl2J21, 



We see from Eq.( |4.72j ) that the two-point of this field decays algebraically, 



c 4 

< q 12 (z,z)q 2 i(0) >oc — 



Z 2 Z 2 



It was recently proposed by Zirnbauer that the PSL(N\N) sigma models describe the 



integer Quantum Hall transition ||16|| . The kind of disorder considered in the present paper, 
which, as we have explained, leads to the PSL(N\N) sigma models, is not generally believed 
to correspond to kind of disorder needed for the Quantum Hall effect, at least not in an 
obvious way. Nevertheless this is an interesting proposal that needs further investigation. 



D. Generalized Random XY Model (Broken Time- Reversal Symmetry) 

In this subsection we discuss a more general random XY model (at the free fermion 
point). This model may be defined in terms of an associated 2-component Dirac hamil- 
tonian of type H 2 , where now all random potentials have both real (A x , A y , V, M) and 
imaginary (A' x , A' y , V, M') parts. As discussed in Appendix B, in this case the corre- 
sponding 4-component hamiltonian H4 lacks in general time reversal symmetry (but still 
has particle-hole symmetry), and the corresponding field theory action possesses GL(1\1; C) 
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global SUSYQ This subsection is devoted to the 1-loop RG equations of this theory. To 
summarize the result, we find that the theory flows off to strong coupling. The strong 
coupling physics should be captured by a non-linear sigma model of the kind discussed in 
Section VI. 

Performing the disorder average over independent random variables m = V + iM and 
H = i(V' + iM') with variance g m = mm* and = fifi* respectively^, we arrive at the 
following lagrangian for the averaged theory 

yJaaJ aa J aa J aa ~ l ~ ^*aa^ 

o=l,2 



^susy ~ 9l E i^aaJ aa J aa^ aa + G aa G aa G\ a G a 

+ 92 E {jbaJab — J'baJ ab + G a bG ab — G\ a Gb^) (4.80) 
a,b=l,2 

where 

9\ = 9m - 9?, 92 = 9m + 9n (4.81) 
The additional contribution from the random vector potentials is 

a^b 

£Ly = 9?(-±) E (J'aa-Jaa)(j'aa-Jaa)+92(-l)J2 (J' aa ~ J aa)(j' bb ~1 bb) (4.82) 
o=l,2 0=1,2 

where 

9i = 9a - 9a>, 92=9A + gA> (4.83) 

[gA and g& are the variances of the imaginary and the real vector potential, respectively] 
We recover the time-reversal symmetric case, Eq.( 4.53j ), by letting g^ = g& = 0. We note 



that in the special case where only the potentials A x , A y , V, M' are non-vanishing, which 
corresponds to g m = 0, g^ > 0, we recover time reversal symmetry [Eq. ( P6|) of Appendix 
B]. The RG equations for this case are identical to those of the case g m > 0, g u = 0, studied 
in the the previous section, upon letting g m — > —g a - This can be seen by replacing in 
Eq.(p8|, |B9] ) the left moving fields ipn by iipii, ipL2 by —iipL2 ( and ijr L1 by —iip L1 , and ijr L2 
by anc l similarly for the /5-fields. All right moving fields remain unchanged. With 

these redefinitions the coupling constants in Eq. (|4.80 ) become g\= g2 = On the other 



hand, as these redefinitions do not change the kinetic term in Eq.( |Bllp of Appendix B, the 



OPE's of these fields remain unchanged. This implies that the OPE's of the gl(2\2) currents 
in Eq. (|4.51|) remain also unchanged, and the entire analysis of sections IV A, B remains 
unchanged upon replacing g m by —g u . 



11 GL(N\N;C) symmetry for the iVth moment averages. 

12 using the action of Eq.(Bll) of Appendix B 
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Let us now return to the most general case. The 1-loop RG equations for the coupling 
constants of the action above are derived in Appendix C ( Eq.'s ( p2|) , ( |C3|) with the result^]: 



<kh_ _ q dg_ 
dl dl 



0, JL- = 2[{g m f + { gil ) 



^ = Dg A ,g 2 , ^- = Ag mgfl (4.84) 



The coupling g 2 flows according to 

d 2 log (g 2 ) „dg A 



dl 2 dl 



D ^7T = D [(92Y ~ (9i 
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which means that g 2 flows to large values. ( Note from Eq. ( [4.8 1| ) that the expression 



K92) 2 — (gi) 2 } is always non-negative.) The strong coupling physics should be described via 
the sigma model technology of section VI. 

V. DENSITY OF STATES OF THE HATSUGAI ET AL. DELOCALIZATION 
TRANSITION FROM CURRENT ALGEBRAS 

In this section we apply some of the above results to the problem of localization of 
electrons hopping with real amplitudes on a square lattice with flux ir per plaquette || . In 
the continuum limit the hamiltonian of this system is the 4-component (hermitian) Dirac 
hamiltonian H4, 

*=Q " 2 ) (5 - 1) 



where H 2 is defined in eq. (|1.4| ). The reality of lattice hopping amplitudes implies that 
the hamiltonian if 4 must have time-reversal symmetry, which is manifest with this form 
of H 2 (see Eq. (p2|) ) . The corresponding SUSY field theory action (Eq.( |Bll| ) of Appendix 
B with n = 0) can be seen to possess GL(2\2; R) symmetry. (This model has also been 



studied numerically in Pb| .) In this subsection we compute the density of states of the 
hamiltonian if 4 . In particular, we will be interested in the eigenstates of the Schrodinger 
equation = E^, where \l/ is a 4-component wave function. The single particle Green 
functions are defined by the functional integral J D^/*D^f exp(— Sq), where 

72 

Si = I ^ ] (x)i (H A - S) 9(x) (5.2) 
and 8 = E + ie. For e = + , the functional integral defines the retarded Green function: 

G R (x, x 1 ; E) ab = lim(x,a|— . \x', b) = lim l i(* a (x)*J(x')> (5.3) 



13 In terms of the notations of the Appendix, D = (d 22 — d^i 
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Making the identifications^ (Appendix B) 

* = V>k V>£1, M*, * f = (V&i> ^Il> fe, 4>L2) (5.4) 

the action takes the form of an £ perturbation of the N = 2 species model: 

72 

S= f — &( x )i (H A - £) (5. 
d 2 x 2 

Y~ S i^LcMLa + ^Radz^Ra + m(x)^ Ra ^ La + m*(x)^[ a ^ Ra ) 
a=l 

+ tylxipLi + V'LV'ia) + A z (Tp R1 ip m + tPm'fpm) ~ iE j (Vii^Io - V'LV'ii + ^ c 

For simplicity, we can set the physical A-couplings to zero from the beginning and then 
let the coupling be generated dynamically under renormalization. Introducing ghosts 
and integrating over disorder leads to the effective action: 

S = Sit 2) (5.6) 

where 

Oe = V44 2 - + PliPm + fiA + 'h.c' terms (5.7) 

(see Appendix B). 

Of interest is the density of states (DOS) 

P (E) = l™ ( tf 4 -E) = ± to ImTY ( 7 ^ 73 ^) < 5 - 8 ) 

where V is the two-dimensional volume. The averaged density of states can be expressed in 
terms of the averaged retarded Green function 

p(E) = - lim Im tr G R (x,x;E) (5.9) 

7T e^0+ 



where tr denotes the trace over matrix indices of the Green's function of Eq.( [5.3| ). 

Thus the density of states is related to the one-point function of the, say, fermionic part 
of the operator Oe which couples to £: 



p(E) = - lim Jte<i(V>LV&i - ^Ri))e (5.10) 

where the one-point function on the right hand side is computed using S e s with the £ term. 

Generally, the £-term spoils the complete solvability of the theory. Our strategy will be 
to use our solution of the £ = theory in conjunction with the renormalization group to learn 
something about the density of states. Let us view Se = —%£ J ^Oe as a perturbation of 



The superscript * denotes the transposed vector. 
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the nearly conformal field theory. The anomalous scaling dimensions of fields are controlled 
by S e ff when £ = 0. 

We need the anomalous dimension T of the field Oe which contains : V'iiV'ja :=: 
exp(-iip L1 + iip R2 ) := exp{-z(0! + 2 )/2} exp{-i(0! - 2 )/2 } H er e, 4> a = Lf La - ip Ra 
denotes the dual field, whose 2-point function is found from Eq.( |4.66| ), (|4.74|) to be 

< fc(z,z)0+(O) > 9A=0 ~ 4g 2 ( 8 log 2 (zz/a 2 ) 

The requirement that for arbitrary g A this function is to satisfy Eq. ( [4.41] ) gives 

< 0+(z,*)0+(O) > 9A ~ 8( 4 g A log(zz/a 2 ) + 4g 2 ( 8 log 2 (zz/a 2 ) 

This together with 



< 0+(z,z)0 6 (0) > 54 ~ -2( 2 5 ab \og(zz/a 2 ) 

t . 

R2 ■ 



yields the anomalous dimension^ of the field : Vxir 



r(g A ) = ±(\l-4( 2 g A ) + h(g) 

where h(g) is independent of g A . Since the action is dimensionless, 

dim(E) = 2-r (5.11) 



ensuring that p(E)dE has dimension 2, as it should in two dimensions. 
The renormalization group equation for p(E, gA) 

( (2 " V)E ^E + Pa{9) Wa ~ r( ^ } ) nE,9A) = ° (5 ' 12) 

implies that the renormalized density of states p(l) = p(E(l), g A (l)) (where e l > 1 is the 
rescaling factor) is related to the unrenormalized DOS p(0) by 

p(0)/p(l) = exp{- / T(l')dl'} = exp{ / [2 - V(l')\dl' - 21} 
Jo Jo 

Here 

dE(l) dg A (l) 

-dT = [2 - T]E > ~dT = ^ 

which permits to express the /-dependence in terms of the the renormalized coupling constant 
E R = E(l) and its bare value E = E(0) < E R . This yields for the DOS 



E 



p(E)/p(E R ) = exp (-2[g A (E R ) - g A (E)]/(3 gA ) (5.13) 



15 For N > 2 the free field property that leads to exact expressions such eq. fl4.43p is not expected 
to hold. However if we assumed that for the purposes of computing the anomalous dimension we 
can parallel what we did for N = 1, we would obtain h(g) = C 2 /2- 
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(valid since /3 gA (g,gA) is independent of gjC). From this we find in the limit E/Er — > 



p(E)/p(E R ) ~ ^ exp{- v ^ (1 + 2g) V log(^)} 

This diverging (but integrable) density of states was not seen in the simulations of 
Rather, a power law varying with the strength of disorder g was found. Such power law 
behavior of the DOS is obtained from the R.G. analysis above for intermediate not asymp- 
totically small energies. 



VI. GADE-WEGNER DELOCALIZATION TRANSITION VIA 

GL(N\N;C)/U(N\N) SIGMA MODEL 

A. GL(1\1;C)/U(1\1) Sigma Model 
Localization problems are usually described in terms of sigma models, using replicas |27 



28fl or supersymmetry p0| . 

In this section we study the Gade-Wegner localization problem [L2]|. A version of this 
model consists of electrons hopping on a 2-dimensional square lattice with arbitrary complex 
hopping amplitudes (and flux 7r per plaquette). The continuum limit of the hamiltonian of 
this model is again, as for the model discussed in section V, of the form of a 4-component 
Dirac hamiltonian. In contrast to the latter model, however, the lattice hamiltonian of the 
present model, and thus also its continuum limit H4 lacks time reversal invariance due to 
the complex hopping amplitudes. As discussed in more detail in Appendix B, the associated 
Dirac hamiltonian H4 has now arbitrary complex random scalar potential ( with real and 
imaginary parts V and V) and complex Dirac mass ( with real and imaginary parts M 
and M') terms, which may be combined into two independent complex random variables m 
and /i, with variance g m and g^ respectively. The explicit form of the hamiltonian is (from 
Eq.([B7|) of Appendix B): 



#4 = H) 



/ 


e 





m — 







e 


d 




— m* + fi* 





€ 


V 


d 


—m — fi 






H d \ 

m* + fi* 


e ) 



(6.1) 



where 



m = (V + iM), fi = t(V' + iM') 



In general, an arbitrary complex vector potential will be generated upon RG flow in the 
corresponding SUSY field theory^. We set the vector potentials to zero initially, and let it 
be generated upon renormalization group transformations later on. In order to study this 
system we use the conventional sigma-model technique. 
The lagrangian is 
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see Eq. (gisj) . 
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where £ = E + ie (e = + ) and 



Let us define new fields 



T 



/ rl 




/ VI 


r2 






^3 








\<Pi) 




/ 


^3 


T+\ 


, T = 


03 


T_ / 


-■02 




\ 


>-02 



Tt = (^, 01, ^, 0I) = (Tt, Tl), 

ft = (^s, -v 2 *, -^) = (n, rl) 

Using Eq. ( |B10| ) of the Appendix the lagrangian can be expressed as 

)f Jt 4 - if (f f T + T f T 



£ = T f £)T - T f £>T + m^aOTjTi 



m± = (m ± /i) 

where a sum over the repeated index i — ±1 is implied and 



£> = ((7+ <g> 1)9, + (cr_ ® l)dg 



( o 



d- z 
\0 



«9 2 





ch 





o / 



When 8 = 0, the lagrangian has the following symmetry: 



T — > £T 



/(G-^t 
I G 



T. 



T -> £T 



G 

(G- X )t 



T 



(6.2) 



(6.3) 



(6.4) 

(6.5) 
(6.6) 

(6.7) 
(6.8) 



(6.9) 



(6.10) 



where G is a 2 x 2 supermatrix with bosonic (fermionic) diagonal (off-diagonal) elements. 
That this is a symmetry follows from Q^{a±®l)Q = G'(a±®l)G = 0±®1 and G'G = 1- Since 
G is only required to be invertible, G is an element of GL(1|1; G), i.e. the complexification 
of GL(1|1). In the presence of the S term, the above transformation continues to be a 
symmetry only if G'G = G'G = 1, which implies G^G = 1. Thus the S term breaks the 
symmetry to i7(l|l). 
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Without loss of generality the random variables (m+, ml) and (m_,ml) may be taken 
to be statistically independent^ with the distribution ( [L.3|) each, and the same variance g. 
Upon disorder averaging one then obtains the term in the effective lagrangian: 

£e// = <?£(rlr 4 )(rjr,) (6.11) 
t=± 

Introduce a grade [a] for the vector index a, where [1] = [3] = 1 (fermions), [2] = [4] = 
(bosons), such that 

T a T 6 = (-l)™T fe T a (6.12) 
Then using the fact that (-1)W 2 + 2 WM = (-l)W, one finds 

(X;f ip )(f *J iq ) = (-l)^(T iq T* p )(f vp %) = StrMM (6.13) 

where Str denotes the supertrace and Mj, are the 2x2 matrices of bilinears: 

M i = T i Tl M i = f i f\; (M l ) pq = T ip T* q , (M,) p? = %f* q , p,q= 1,2 (6.14) 

Let us introduce two 4x4 block diagonal supermatrices of fields Q, Q 

n= (Q+ o \ h= (Q+ 

v_ I o qJ' ^" v o g_ 



and define 



M_ / ' V M_ 



Consider the following lagrangian: 

C = Str ((£> - gQ - i£)M - (D + gQ + i£)M - ^QQ) (6.15) 

A functional integral over the Q±,Q± fields gives the effective interaction £ e // = 
gStr (M + M + ) +gStr(M_M_). If instead we first perform the gaussian functional integral 
over the T, T fields, we obtain the effective action for the Q, Q fields: 

S Q £ = STr log (p- gQ- iS) + STr log {D + gQ + i£)-g J StrQQ (6.16) 

In the above equation, STr incorporates the integral d 2 x: For a diagonal functional operator 
A(x, y) = A(x)8(x — y) we define STr^4 = 4y / ^StrA(x), where a is an ultraviolet cutoff. 

The symmetry transformations on M, M follow from (|6.10|) . The lagrangian (|6.15| ) then 
has the following symmetry when 8 = 0: 



17 For general values of g m and g^, the variables m + and m_ will be correlated and will have 
different variances. Changing the ratio g m / ' will however not affect the universality class of the 
sigma model. Therefore we may chose this ratio to be unity. 
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M^QMQ\ M^QMQ ] (6.17) 



When £ 0, the symmetry is as in ( 6.17 ) with G an element of C/(l|l). 

Due to the linear terms in S Q 5, the fields Q, Q develop vacuum expectation values (Q), 

(Q). Let us take the vacuum expectation value (VEV) to be ((Q±) a b) = q±8 a b, ({Q±)ab) = 
q±S a b- Minimizing Sqq with respect to Q,Q, one finds that for £ — > 0, the VEV's q = q± 
and q = q± must be solutions to a set of self-consistent equations which have the solution 

q = q= — (exp(— )-lY l/2 (6.18) 
a 9 9 

which for small g goes to zero exponentially as follows: 

1 2n 

g = g~ — exp( ) (6.19) 

a 9 9 

This follows from the saddle point equations 

Q = - <x\(D-gQ)- 1 \x>, Q = < x\(D + gQ)~ x \x > 
which, when written in momentum space, lead to 

£f l0g(l + 1 ) = g 

4tt ( a gqf 

(and a similar equation with q <->• q). 

Since the VEV breaks the GL(1\1; C) symmetry to Z7(l|l), there are massless Goldstone 
modes. These Goldstone modes as usual can be viewed as a symmetry transformation of 
the VEV. Thus, let us define the fields 

■U-i,n\r-i - fl GG * \-(lU 



« = (6T'<Q)6- = ^ ?(GG , r iJ = ^ qU -i) (6-20) 

The latter are expressed in terms of U = GG\ An element G G GL(1|1; C) can be factorized 
as G = G a G u where Glfiu = 1 and G\ = G a . One has U = G a G u G\jG a = G 2 a1 thus the field 
U lives on the coset space GL(1\1;C)/U(1\1). 

A parametrization of the Goldstone modes U which manifests the required properties 
that U be invertible and U' — U is the following: 

U ~{0 e^l V2x f 1-xxVVO e*») ^ Zi) 
rj-x = ( e-^ \ (I - X ] X -V2 X \ ( e-^ 
V e-*VV— N/2x f 1-XXVV e"* 

The fields </>i 2 are bosonic whereas \i X are fermionic. 
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The low energy cr-model for the Goldstone modes U is obtained by substituting Q,Q — > 
U, U and performing a derivative expansion using STr log (A + D) — STr log A + STr log(l + 
A^D). The zero-th order term in £ is 



2aV 



Str ( DU^DU- 1 + DU- l DU~ l 



q 



a 9 



(6.22) 



Keeping only the linear term in £, and expressing everything in terms of U, one finally 
obtains 



S 



d 2 x 



1 



8ttA 2 
where for small g 
1 

A 2 ~ ; " 



Str (dM-'dM 



(Strf/- 1 ^^) 2 - i£\ 3 Str([/ + E7- 1 ; 



a 2 g 2 q 2 



Ait 

4exp( — ) > 1, 
5 



2 2 . 2ti s 
— — w - exp(— 
«S9 a g 



As explained below, the A^/A 2 term is generated under renormalization. 
In terms of the parametrization ( |6.21| ) we find the following expressions: 



(6.23) 



(6.24) 



str (djj-%u) HW^ + ^Al + A 



- (StTU~%U 



$tr(U + U- 1 )=2(1- xx f ) cosh 2<p 2 - 2(1 - X j x) cosh 2<^ 



(6.25) 

(6.26) 
(6.27) 



where <j> ± = (fa ± 0i) (6.28) 



(The effect of a Wess-Zumino topological term is discussed below.) 

We may study the sigma model as a perturbed conformal field theory. As an example we 
will compute the density of states for the present model. Rescaling the fields 0i 2 — > j<t>i,2, 
X — *■ t Xi the sigma model action can be written as 



r d 2 x f\ 2 \ 
S = S CFT + j — i-O g + X A 2 - i£0 E \ (6.29) 

where the conformal field theory has the same field content as our <7/(l|l) model, with action 

d 2 x 



and 



S cft = J ((d.fo) 2 - (d^) 2 + d, X ] d,x) (6.30) 



O g = \x ] X (d^-) 2 (6.31) 
2 = - {d^-f (6.32) 



O e = 2A 3 



A 2 A 2 
[1 - -^XX ] ) cosh(A0 2 ) - (1 - -jX ] x) cosh(A0i) 



(6.33) 
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Since x X is a logarithmic operator p9[, it mixes under RG scale transformations with 



the identity operator, thereby generating the A a coupling. 

This can also be seen, alternatively, by making the following change of variables 

which yields 

- -Str (d,U- l d,u) ={8^+8^ + e- 2 *- 8^8^} - 8,[e^- (6.34) 

The last term can be eliminated by shifting </> + , 

+ - $ + = + - e~ 2 *- (6.35) 

Let us now also consider a topological Wess-Zumino-Witten term with coupling k. To lowest 
(cubic) order this term is easily computed with the result 

S wzw = 4- / d 3 x e ijk Stilu^diU [U- l 8jU, U~ l 8 k U]X = (6.36) 
2An Jm 3 I J 

% k f 

= — d 2 x (du4>-)e^ v [{8 u £})(i — ^(<9„£)1 + quartic and higher terms 

Z7T J 

where, as usual, M 3 is a three-dimensional manifold whose boundary is the two-dimensional 
space of interest and e^ u the two-dimensional antisymmetric tensor. 

The sigma model including the WZW term can be solved using the Lagrange multiplier 
method of section IV. In order to see this, note that the fundamental matrix field is of the 
form 

U = e*+ V (6.37) 



where V does not depend on $ + . This implies that the commutator in Eq.( |6.36D does not 



contain any $ + -dependence. Since the supertrace of a commutator vanishes, the WZW 
term does not contain any <3> + dependence either^. Therefore, as for the current-current 
perturbation of section IV, $ + acts as a Lagrange multiplier, and we can use the same steps 
to solve the present theory. 

Making the same rescalings of the fields as above, the action for the sigma model with 
WZW term is 

S = Soft + J ^ Q O g - h\ 3 O wzw + X A 2 - iEO^j + ... (6.38) 

where^ 



18 Note that this argument is equally valid for the GL(N\N)/U(N\N) generalization of section 
VI B below. 

19 No confusion should arise between the fields <&±, $i 5 2 defined here and those of Eq. (^2,3.3) 
and in Appendix A,B. 
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$ + = ($ 2 + $ 1 ) ; $_ = 0_ = ($ 2 - $ x ) (6.39) 

and only terms of linear order in <3>_ have been writtenF]. Here Scft is given as before by 
Eq.( |6.30| ) upon replacing, ± — > <3> ± , x — * £> an d 

^ = -^-(5^0 (6-40) 

Owzw = \*- i^id^dvt) (6.41) 

2 = -(^$_) 2 (6.42) 

cosh(A$ 2 ) - cosh(A$ 1 ) - (e A$1 + e A * 2 ) 



= 2A 3 



(6.43) 



We proceed with the 1-loop R.G. equations. One has the OPE's: 



O g (z,z)O 9 (0)-^O 2 (0) + -- (6-44) 

O wzw (z,z)O wzw (0) ~ £^ 2 (0) + .... (6.45) 
Thus, as in section IV, to lowest order the /3eta-functions are 

(3 X = 0, A = 0, p x =?L[i-(k\ 2 ) 2 ] + ... (6.46) 

16 

The sigma model coupling constant A, and the WZW coupling k generate A^ terms (denoted 
qa in section IV) of opposite signs. The 1-loop R.G. equations have the same structure as 
those found for the g/(l|l) current- current perturbation model in section IV, where 2 is 
analogous to Oa- In the presence of a WZW term, the sigma model becomes the conformal 
g/(l|l) WZW model at a particular finite value of the sigma model coupling 1/A 2 = k. 

We may solve the theory in Eq.( |6.38| ) exactly using the Lagrange multiplier method of 
section IV A. In order to compute correlators of $ + one needs 

= - (« 2 )]3rf*£« + [1 + (*A 2 )]£«fcU 

As in section IV A this gives the two point function of $ + fields: 

A 2 1 

< d g dg$+(z, z) d w duj$+(w,w) >\ A =o= — - (^A 2 ) 2 ] ■—— — — - + 5 - functions 

2 (z — wY(z — wr 



20 As in section IV, the term linear in 3>_ will be sufficient in the following. 
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Furthermore, as in Eq. (|4.35 ) we find 

< $+(z,z) $_(0) >= -21og(*z), < $-(z,z) $-(0) >= 
As in section IV A, this yields the correlation functions 

< ^(z, z)^(0) > = (1 - 4X A ) \og(zz/a 2 ) - ^[1 - (kX 2 ) 2 ] \og 2 (zz/a 2 ) (6.47) 

lb 

< ® 2 (z, z)$ 2 (0) > = -(1 + AX A ) \og(zz/a 2 ) - ^[1 - (&A 2 ) 2 ] log 2 (^/a 2 ) (6.48) 

16 

as well as the exact beta function (identical to the 1-loop result): 

!£_ ft, _£[!_(«*)•] (6.49) 

We now proceed with the calculation of the density of states. The most relevant operator 
in Oe is cosh(A$2), whose anomalous dimension is 

r(A A ) = dim (cosh(A$ 2 )) = -A 2 (l + AX A ). (6.50) 

The analysis of section V applies to the density of states with X A playing the role of g A , 
yielding 

E f 



p{E)/p{E R ) = -f exp (-2[X A (E R ) - X a (E)]/(3 Xa ) (6.51) 

(valid since f3\ A is independent of X A ). In the limit E — > this gives a divergent (but 
integrable) density of states 

1 4\/2 / 

~ P(E) ~g exp{ - Ayi-(^)> e(£ * /£>} 



in agreement with the perturbative result obtained by Gade |T2| using replicas when k = 0, 
i.e. in the absence of the WZW term. 

We end this section by commenting on the connection with the scale invariant PSL(1\1) 
sigma model. Eliminating the two bosonic coordinates </>± the fundamental PSL(1\1) sigma 
model field becomes 

V ~[ y/fy (l + xh)J (6 - 52) 

whose dynamics is governed by the quadratic action obtained from Soft of Eq. ( |6.30| ) by 
eliminating 0i,02- I n subsection C we briefly discuss the generalization to the less trivial 
PSL(N\N) case. 
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B. GL(1|1)/U(1|1) Sigma Model as a Current-Current Perturbation 

Upon comparing the correlation functions of the sigma model (Eq.( |6.47D ) with those of 
the N = 1-species XY model (Eq.( |4.39| , |4.40|) ), which is a current-current perturbation, one 



finds that the two models have identical correlation functions if one identifies the coupling 
constants as follows: 

2X A = C 2 9a (6.53) 
h 2 [l - (fcA 2 ) 2 ] = (2s) 2 C 6 , (where C 2 = 1/(1 + 2g)) (6.54) 

The identity of the sigma model with the current- current perturbation of the conformal 
current algebra is likely a special case of a more general statement (see next subsection). 

C. GL(N\N)/U(N\N) Sigma Model 

The generalization of the sigma model of the previous section to GL(N\N) is immediate. 
We start with a theory of 2N bosons and 2N fermions (in the previous section N = 1). The 
construction of the previous section provides us in the standard way with a field variable 
on a manifold that we denote by U G GL(N\N;C)/U(N\N). Denoting the elements of 
the Super-Lie algebra gl(N\N) by J± and J a, a parametrization of the fundamental matrix 
field, say of the formQ 

U[X ± ,X A ] = e {x+J+} exp{X_J_ + £ X A J A } 



yields an action generalizing Eq. (|6.25 ). The PSL(N\N) version is obtained by setting the 



coordinates X± to zero. This sigma model in the presence of a WZW term with coupling k 
becomes the gl(N\N) WZW conformal field theory at a finite value A 2 = 1/k of the sigma 
model coupling constant. The latter is a gl(N\N) current algebra at level k. Generalizing 
the observation made in the preceding subsection for the N = 1 case, we suggest that the 
line of fixed points arising from the current-current perturbation of the PSL(N\N) variant 
of the current algebra analyzed in section IV is an alternative and 'dual' description of 
the PSL(N\N) sigma model with WZW term, as the sigma model coupling is varied away 
from the scale invariant point. We plan on reporting in more detail on this connection in a 
subsequent publication. 



VII. CONCLUSIONS 

In this paper we have studied current-current perturbations of gl(N\N) super-current 
algebras. The existence of two quadratic Casimir invariants leads to a 'nearly conformal' 



21 



This follows as in the N = 1 case from the decomposition G = G a G u , see below Eq. (|5.20 ). 
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theory for any value of the coupling. Current algebra techniques allowed us to compute 
exactly the correlation functions of all current operators in the perturbed theory, as well as 
the exact /3etaf unctions (in a particular RG scheme). The current correlators are conformal 
except for correlation functions involving the currents corresponding to the trace and the 
supertrace of the super-Lie algebra. In the case N = 1 we have computed the correlation 
functions of all fields. We have applied this to two different, but related disordered models, 
the 2D random field XY Statistical Mechanics model, and the derealization transition of 
electrons hopping on a 2D square lattice with 7r-flux per plaquette and real hopping ampli- 
tudes (HWK, ||). Both theories can be formulated in terms of the 2D Dirac hamiltonian 
subject to random real mass, imaginary scalar and imaginary vector potentials. This hamil- 
tonian is invariant under two discrete symmetries, charge-conjugation (particle- hole) and 
time-reversal, in every realization of disorder. For the random XY model we have com- 
puted all correlation functions involving lst-moment averages (described by the '1-species 
theory'), as well as certain correlation functions involving Nth-moment averages (namely 
those involving the gl(N\N) currents). In particular, we have computed the correlation 
function of the 'Edwards- Anderson order parameter', which is scale invariant and of scaling 
dimension x = 2, in the 2nd moment (N = 2-species) theory. For the random hopping 
model of Hatsugai et al. we have obtained the density of states which shows a divergence at 
zero energy. In section IV D we derived the 1-loop RG equations for a generalized random 
XY model: When formulated in terms of a random 2D Dirac Hamiltonian, the latter still 
exhibits particle-hole symmetry, but lacks time-reversal symmetry, as opposed to that of 
HWK. The derealization transition exhibited by this more general random Dirac hamil- 
tonian is in the Gade-Wegner universality class, describing hopping of electrons on a 2D 
square lattice with general complex hopping amplitudes. In section VI we have derived a 
SUSY sigma model from the underlying random Dirac hamiltonian, found to be defined on 
a target manifold which we denote by GL(N\N; C)/U(N\N). We have solved this sigma 
model exactly for N — 1 (including also a WZW term), and shown that it is identical to a 
current-current perturbation of a GL(1|1) super-Lie current algebra. The density of states 
is obtained from the sigma model. 
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APPENDIX A: SUSY APPROACH TO THE RANDOM FIELD XY-MODEL 

In this Appendix we show that the SUSY approach to the random field XY model leads 
to the N = 1 species action discussed in the bulk of the paper. 

The non-hermitian hamiltonian relevant for the 2D random field XY model is 

H 2 = {-id x + iA x )a x + {-id y + iA y )a 2 + Ma z - iVl 2 
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In order to be able to integrate over bosonic variables, we need to consider the hermitian 
hamiltonian 



Ha = 







H 2 




(Al) 



= -%d x o x ® n - id y a 2 ®n + A x a 1 ® r 2 + A y o 2 ®r 2 + Ma 3 ® n - Vl 2 <8> r 2 (A2) 

In the case of interest where all potentials A x ,A y , M, V are real, Ha is time reversal invariant, 
which means that H± can be brought into real, symmetric form. Explicitly, this can be 
achieved by conjugating H± with the matrix 



where 



L7 CT = (1 - ia 1 )/V2, U T = (1- ir 3 



yielding the real symmetric hamiltonian 

H S A = U^H 4 U = id t (a! <g> r 2 ) + id y (a 3 <g> r 2 ) - 7(1 <g> r x ) - M(a 2 ® r 2 ) (A3) 

For integration over fermionic variables an antisymmetric form is needed. This can be 
obtained by conjugating by 

U' = U a ® 1 2 



HI = U'^HaU' = vdjpx ® n) + id y {a z ® n) + V(l ® r 2 ) - M(a 2 ® n) 



(A4) 



All correlation functions relevant for the XY model can be obtained from an action 
constructed with a (real) 4-component fermion field x, and a real 4-component boson field 
(p, defined as follows: 



X 



fXi\ 

X2 
X3 



with Lagrangian (density) 
C = Cf + Cb-, 



Xa , 



Cf^tfUtHtU'x, 



(<Pi\ 



( v+ 
\<p- 



It will prove convenient to introduce 4-component complex fields and their adjoints by 



Fermions : 
Bosons : 

so thatQ 
Cf = iWHaV = 2i^ + H 2 y 



* = 



X 



$1 



C h = i&H d <5> = 2i$t# 2 $_ 



2mlHt<$>, 



22 Note that the latter equalities involving ty±, &± arise because i/| is symmetric, and because 
V and U do not mix the 2x2 blocks involving H 2 . 
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The so-defined fields are not independent from their adjoints. Rather, the above definitions 
immediately give the following relations: 

(F)* = u'*x = U'*U*V, ($ t )* = <F = U*rf$ 

Since U'*U'^ = U*W a = = and U*W = -a l <g> r 3 this implies 



^3 
\1>J 



in 







(-<P*\ 


















\ ^3 / 



Therefore we may chose 
as well as 



and 



and 



^3 
^4 



Pr 



(AS) 



(A6) 



as our independent 4 fermionic and 4 (real) bosonic integration variables. The resulting 
super symmetric action 



£susy=*(^ r , 4>[)H 2 (^ L R )+i(Pl /3l)H 2 



Pl 
Pr 



is the N = 1-species theory discussed in the bulk of the paper. 



APPENDIX B: THE MOST GENERAL PARTICLE-HOLE SYMMETRIC DIRAC 

HAMILTONIAN 

In this appendix we discuss in some detail the SUSY invariant action for the most 
general random particle-hole symmetric, but not necessarily time-reversal symmetric Dirac 
Hamiltonian. We start from the hermitian quantum mechanical Dirac Hamiltonian iJ 4 
defined by 

where 

H 2 = (-id x + iA x )o x + (-idy + iA y )a 2 + Ma 3 - iVl 2 

(As usual, infinitesimals e — > + are needed to extract quantum mechanical Green's func- 
tions). This model is time reversal invariant if all the potentials (A x , A y ,V, M) are real. In 
this case the time reversal operation is implemented by 

THIT = H 4 , with : T = a ± <g> r 3 (B2) 
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We may consider a more general (non-hermitian) model, where we add to any of the 
potentials an imaginary part, A' x ,A' y , M', V . Note that A' x , A' y corresponds to a real vector 
potential. The most general hermitian hamiltonian Ha incorporating all these potentials is 

H 4 = -id x a 1 <g> Ti - id y a 2 ® n + A x a 1 <g> r 2 + A y a 2 <8> r 2 + Ma 3 r x - Vl 2 ® r 2 (B3) 
- A' x ai ® ri - A' y er 2 ® ri - M'cr 3 ® r 2 + V'l 2 ® ri (B4) 

In general, this hamiltonian is no longer time reversal invariant^. However, it is always 
invariant under particle-hole transformations 

C HiC = -Hi, with C = 1 ® r 3 (B5) 

Time reversal symmetry is recovered for the special case where only A x ,A y ,V, M' are non- 
vanishing. In this case the time reversal operation is 

T'HlT' = H Al with: T = a 2 ®l 2 (B6) 

These time reversal symmetry properties are summarized in the following table: 



non-vanishing potentials 


A x ,A y ,V,M 


A x ,A y ,V',M' 


others 


time reversal 


YES (T) 


YES (T') 


NO 



Let us now return to the most general situation, where time reversal symmetry is in 
general absent. The 4x4 matrix corresponding to the hamiltonian Hi has the following 
block structure: 



% H A 



where 
and 



/ e m-fi d + A z + iA' z \ 

e d + A* z + iA'* z m* + 

-m* + /j* d - A z + iA' g e 

\d-A* z + iA'* z -m-fi e / 



m = (V + iM), 



(-fi) = (M f - iV) 



A z = (A x -zA y ), A Z = A* Z (-i) A ' z = (A' y -iA' x ), A'- = A' 

This matrix is manifestly anti-hermitian as it shouldFJ 



(B7) 



23 the imaginary parts A' x ,A' y ,M',V' break invariance under time reversal, using T defined 
above. 
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note that d cx [d x — id y ], d cx [d x + id y ], so that = [— d x — id y ] = —d. 
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Averaged Green's functions of the hamiltonian H4 can be obtained from the following 
action. We introduce^ a 4-component complex bosonic field $ and its adjoint $' , 



$ = 



(<h\ (Ph\ 

Pr2 
Pli 



\PbxJ 

= 02, <S <f>t) = (P]tL, Pll, -PR2, 

as well as a 4-component (complex) fermionic field \I/ and W: 



(B8) 







M 2 \ 






%2 


^3 






W 4 y 







* f = (v>i, vs, = (A, v>L Vw) 

The SUSY lagrangian (density) reads, in the absence of the vector potentials: 



(B9) 



+ 



V^WVWi) + (0t<90 4 +01<90 1 

171+ f^3+0I03Vm-te(-^2)+01(-02 



+ e 



^(-^)+(-^)^ 3 j+( v 035(-0 2 )+(-0;)^ 3 

^^1+0301) +m- ((-V'2*)V'4+(-02)04 
V^l + 0101 + (-V2) (-^2) + (-0D("0 2 ) + «3 + 0^03 + ^4>4 + 010J (BIO) 



+ 



+ 



where 

m + = (m + fi), m_ = (m — /x). 

This may also be expressed in terms of the R- and L-moving fields, defined above: 

r = 

^susy 

= £(4a^a + /4^a + Aa^La + P{ a dP L a) + 
a 

+/" ((V'mV'Li + 4i^0 - (Afe + PkiPl*))-^ ({Ai4>Ri + PliPri) - (Ai^Ri + PmPRi) 



+ 



25 The symbols /3k, /?|^ will be used to exhibit the chiral (L/R) nature of the corre- 

sponding fields (see below). 
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+ e (^m^L2 ~ V>LV>li + ^L2^m - ^Life + P R1 p{ 2 + /4a/4i - #m&2 + PmPn) (Bll) 
The vector potentials give rise to an additional contribution to the action of the form: 

lA' z (V4lV>Rl + PhPm ~ TPR2*PR2 ~ PhPxt) + lA'z (ipljIpLi + P j L1 P L l ~ ^L2 ~ P1 2 Pm) = 

(J 11 - Jn) + (J'22 ~ J22)] + AgftJ'n - J n ) + (J' 22 - J 22 )] + 
J '11 - Jn) - (J'22 - J22)] + iA'z[(J'u - J u ) - (J'22 - J22)] (B12) 



+ 
= A Z 

+ %A' 



APPENDIX C: 1-LOOP RG EQUATIONS FOR THE GENERALIZED RANDOM 
XY MODEL (BROKEN TIME REVERSAL SYMMETRY) 

In this Appendix we analyze the 1-loop RG equations for the random XY model without 
time reversal symmetry. 

Let us first consider the RG equations for the couplings gi = g m — <? M and g 2 = g m + g M , 
defined in Eq.'s ( |4.80| , 4.81 ). The most general formP] of the 1-loop RG is: 



-Jr = a i9i + b i9i92 + c x g 2 2 + (d\ x g x + d\ x g 2 )g x + [d\ 2 g x + d\ 2 g 2 )g^ 

d-92 „ 9 1 1 1 „ „2 , / t2 „ , t2 „ \A , / t2 „ , t2 „ \ 4 



We proceed in the following steps, to simplify these equations: 

Step 1: We know that both RG equations vanish when g 2 = g 2 = describing two decoupled 
1-species theories. This gives a\ = d\ x — a 2 — d\ x — 0. 

Step 2: Adding and subtracting the two equations: 

d{91 * = = (fol+62)^l^+(ci+C 2 )^+(4 1+ ^jW+[(^2+^2)^l + (42+42)^]^ 



d{9l dl 92) = (-2)^ = (6 1 -6 2 )^^+(c 1 -c 2 )^+(4 1 -4J^^+[(^ 2 -^)^+(42-4 2 )^]^ 



26 Note that we have omitted terms on the r.h.s. which are quadratic in the vector potential 
couplings, i.e. of the form gfgf. The operators coupling to gf are of the form J aa and J'bb (see 



Eq.( 4.82D ). Since there are no poles in the OPE of such operators with each other (see Eq.( |4.5lD ) 



they do not generate 1-loop RG flows. 
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Setting = g& = 0, both equations must vanish identically. In this case we have 
9i = 92 = 9m and 9\ = 9% = 9a-, which yields, when inserted into the RG equations: 

h + ci = + d\ 2 + d\ 2 = b 2 + c 2 = d\\ + d\ 2 + d\ 2 = 

Step 3: Finally, setting g m = g& = 0, again both equations must vanish. In this case we have 
gi = —g^g-i = g^ and gf = g 2 — g A , which yields, when inserted into the RG equations: 

—bi + ci = d\ x - d\ 2 + d\ 2 = -b 2 + c 2 = d\ x - d\ 2 + dj 2 = 



Combining Steps 1, 2 and 3 we find: 



a i 



bi 



= b 2 ~- 

(4l 



Cl 



d\i 

= c 2 



a 2 
0, 



d 2 n = 
d\ 



d\ 2 



d\ 2 ) 



(d 2 21 + d\ 2 ) = 



(Cl) 



The 1-loop RG equation thus simplifies to 



^ = + d l 22 gf)g 2 



dl 

Using Eq. (|4.83|) this becomes 



-Or = (d 2 21 gi + d\ 2 g 2 )g 2 



dgi 
dl 

dg 2 
dl 



-(4i - d\ 2 )g A ,g 2 



-{d\ x - d 2 22 )g A >92 



Let us now consider explicitly the relevant OPE's: The coupling constant qa> couples to a 
left-right bilinear of a Cartan generator. Since the simple pole term in the OPE of a Cartan 
generator with any operator can only give back the same operator, times a number (= the 
eigenvalue of the Cartan generator when applied to this operator), we know that 9a>92 can 
only generate again the operator coupling to g 2 . Therefore we know that (d^i — ^22) = 
and the RG equations simplifies further to: 

^ = -(4 - d 2 22 )g A/ g 2 (C2) 
(Recall that g 2 = g m + g M , from Eq. (|4.81| ).) 
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The corresponding RG equations for the vector potential couplings are of the form 

dg A 



dl 
dg A , 



<*ig m + PiQmgn + 7i0 P 
+ fang? + I19I 



dl 

Since g A > is not generated when either g^ = or when g m = 0, we have a 2 = 72 = 0. In 
addition, the flow for g A is the same in both cases, implying ai = 71. Adding and subtracting 
the resulting equations 

dgf d(g A - g A i) . 2 2 . 

~dT = — — = ai ^ m ^ ^ ~ P 2 > 9m9 » 

dg£ d{g A + g A <) 2 2 . 

= ^ = aiUm + V + (#1 + P2)g m g^ 

When (yf 2 = (i.e.: = —g m , 9i = 2g m ) we know that 

dgt 2 dg£ 

This finally yields the following RG equations (fii = 0, (3 2 = 2ai): 

^ = 2[(, m ) 2 + (^) 2 ] ; ^ = 4 ^ ( C3 ) 

where we have used that «i = 2 [see Eq.( |4.1UD ]. 
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